Spectrum of an electron spin coupled to an unpolarized bath of nuclear spins 
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The main source of decoherence for an electron spin confined to a quantum dot is the hyperfine 
interaction with nuclear spins. To analyze this process theoretically we diagonalize the central spin 
Hamiltonian in the high magnetic B-field limit. Then we project the eigenstates onto an unpolarized 
state of the nuclear bath and find that the resulting density of states has Gaussian tails. The level 
spacing of the nuclear sublevels is exponentially small in the middle of each of the two electron 
Zeeman levels but increases super-exponentially away from the center. This suggests to select 
states from the wings of the distribution when the system is projected on a single eigenstate by a 
measurement to reduce the noise of the nuclear spin bath. This theory is valid when the external 
magnetic field is larger than a typical Overhauser field at high nuclear spin temperature. 
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Spin dynamics in semiconductor nanostructures has re- 
cently become a topic of great interest due to the possibil- 
ity of using the spin degree of freedom instead of charge 
in electronic circuits [jy] and equally important due to the 
proposal of using electron spin in a semiconductor quan- 
tum dot as a fundamental building block of the quantum 
computing device [21] • GaAs quantum dots are the main 
candidates in practical realizations of these proposals due 
to the well developed manufacturing technology. How- 
ever, unavoidable inhomogeneous hyperfine interaction 
of electron spin with many nuclear spins of the host crys- 
tal acts as a noisy environment that is the main source of 
dephasing for the electron spin at low temperature when 
relaxation due to the phonons is ineffective. 

The limit of fully polarized nuclear spin bath was an- 
alyzed exactly in [3|, including spectral properties. How- 
ever, it is rather hard to achieve a significant polarization 
dynamically, and thermodynamic polarization, requiring 
sub-milli Kelvin temperatures [4] , is still out of reach for 
semiconductors. Currently, a more promising route is 
to actively reduce the distribution width of the nuclear 
Overhauser field by projective measurements d-Q. This 
has been partially achieved in experiments leading to sig- 
nificantly longer decoherence times To further op- 
timize projective measurement techniques it is essential 
to gain a better understanding of the spectral properties 
of the unpolarized system which, so far, have only been 
understood qualitatively. 

In this paper we diagonalize the central spin Hamil- 
tonian for a quantum dot in the high magnetic B-field 
limit using a l/i?-expansion. Projecting the eigenstates 
on an unpolarized state of the nuclear spin bath we find 
that their density has Gaussian tails. Correspondingly 
the level spacing of the nuclear spin sublevels, which is 
exponentially small with the radius of the quantum dot 
in the middle of the two electron Zeeman levels, becomes 
super-exponentially large with detuning away from the 
center, see Fig. [T] This suggests using a finite detuning 
from the bare electron Zeeman energy when one elimi- 
nates the effect of the nuclei by the projective measure- 
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Figure 1: Numerical evaluation of v {E) using Eq. ([2]) on a 
course scale - thick line and Eq. ([5| - thin line {Sj = 1, ro = 8, 
A'' = 18, a fixed external B-field), ^ is a maximum Overhauser 
field, Eo is a shift from Eq. (5). Insets show v {E) on a fine 
scale in the middle of the upper electron Zeeman line and at 
a finite detuning, the average level spacing d was evaluated 
using Eq. ((6)1. 



ment technique 

Our theory is applicable when the external magnetic 
field B is larger than a typical Overhauser field at high 
nuclear spin temperature due to fiuctuations -Bfiuc — 



Ay S/N/jj,, where A/^ is the maximum Overhauser 

field, N is the number of nuclei under the electron en- 
velope wave function, and 5* is a number of degenerate 
hyperfine couplings. At low field B < Bfluc the spectrum 
can be obtained by a numerical solution of the Richard- 
son equations 01 where the 1/B -expansion of the present 
paper can be used as a benchmark for complex numerical 
procedures. 

The spin of an electron in a quantum dot couples to 
nuclear spins in the presence of an external _B-field as 
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H = fiBS^ 



J2A,So 
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where /i = g^B is the electron magneton (in the follow- 
ing we neglect the nuclear Zeeman splitting) , , Sq ~ 
Sq ± iSg are electron spin-1/2 operators and Sj {j > 1) 
are spin operators of nuclear shell j with the maximum 
angular momentum Sj > 1/2 constructed out of 25'^ nu- 
clei of spin-1/2 which have the same hyperfine coupling 
to the electron spin, Sj = X]i=i Ij* i where i labels indi- 
vidual nuclei within the shell, Iji are nuclear spin-1/2 op- 
erators, and N is the number of nuclear shells. Assuming 
harmonic confinement of the electron in all spatial direc- 
tions the couplings are Aj = Aq exp {—rj /r-g) , where Aq 
is the coupling in the middle of the quantum dot, and rg 
and rj are spatial size of the quantum dot and radius of 
jth shell in units of the lattice parameter. 

In ID only two nuclei have the same coupling ignoring 
the isotope effects and assuming equidistant lattice sites 
rj = j, thus the maximum total angular momentum is 
Sj = S — 1. In 2D degeneracy of the couplings gives 
Sj — S — A but the radii of the sequential shells are not 
equidistant because the number of nuclei grows linearly 
away from the center. We thus model the system as a 
set of concentric nuclear shells, r-,- = r + 4m/ (nr) and 
also change the summation indices in Eq. (1), J27=i 



N,Trr/i 



|12| . In 3D the degeneracy is larger than in 



2D, Sj = S = 12, and the number of the nuclei grows 
quadratically away from the center, rj = r + 6m/ (7rr^), 
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This model conserves the number of excitations 
[H,Jz] = 0, where = ^jLo^j^ ^^^'^ total angu- 
lar momentum of each nuclear shell [H, S|] = 0. All of 
them also commute with each other, [J2,S|] = and 
[Sf,S^] = 0. Thus the Hilbert space is partitioned into 
a set of disconnected subspaces labeled by the follo wing 
quantum numbers: n is an eigenvalue of Jz and Ij [l3| 
correspond to S^, S| l^") = lj{jj 1) |*). The latter be- 
comes trivial when all of the nuclear spins have different 
couplings as for spin-1/2 operators = 3/4 is a number 
but is nontrivial when Sj > 1/2. 

The diagonalization in each subspace can be performed 
using degenerate perturbation theory when the B-field is 
large. Splitting the Hamiltonian into the unperturbed 
part Hq = (J^BSq and a perturbation V — "^jAjSo ■ 
Sj defines two electron Zeeman levels, E = but 
leaves the nuclear spin sublevels hugely degenerate in the 
zeroth-order approximation. The latter degeneracy has 
to be lifted via a diagonalization of the perturbation V . 

In the basis of eigenstates of Jz, \'^) = \±,{lj,kj}), 
(^l^f) = 1, y is a diagonal matrix within both of the 
electron spin subspaces where the spin-flip part of V that 
couples opposite electron levels can be neglected when 
the external field is very large. Here ± refers to the "up" 
and "down" electron Zeeman levels and kj are the num- 
bers of nuclear spin excitations on each shell such that 
the quantum number n = (1 ± 1) /2 -|- X^^Li The sec- 
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ond order correction to the eigenenergies are due to the 
spin-flip part of V. Using the matrix elements of V in 
the basis of eigenstates of Jz we obtain 
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i=i 



A. 



A] {21, 



2 , 



(2) 



where the energy denominator in the last term was also 
expanded up to the leading order in l/iiB. Including the 
first order corrections to the eigenfunctions we get 



A A 



\^) = |±, {i„k,}) ± }_J It, {ij,k,}) . 

m=l ^ 



(3) 



The large magnetic field expansion has different condi- 
tions of applicability for the eigenenergies Eq. Q and 
the eigenstates Eq. in the subspaces of unpolar- 



ized nuclear spins kj 



The subleading terms in 



Eq. ^ are small in all subspaces when B ^ -Bfluc 

where Sfluc — \J^^=i ^"j^j //^- But the next (second) 
subleading correction to Eq. is small only when 

B ^ -Bmax where iJmax = Jin^o/'^t^ ill ID and 2D 
(Smax — rQAo/^/Seii in 3D) [1^ is a much larger field 
than -Bfluc- The latter signals that the choice of the 
eigenfunctions, \'^) — \zL,{lj,kj}), is a poor zeroth 
order approximation in the intermediate field regime, 
Bfluc ^ B <C Bmax- The correct approximation can be 
identified by merging the inner nuclear shells with dif- 
ferent couplings up to the radius f = rg-^/ln (-Bmax/^) 
(in units of the lattice parameter) in ID and 2D (f = 

ro (l + Vlii(^max/B)) /V2 in 3D) ^ into a single 
shell with the same coupling Aq. Then, diagonalizing 
H + V, where V = Ej:.,<f (^i " ^i) So • S,, when 
-Sfluc B -Binax instead of the original model H we 
obtain the same result as in Eqs. ^ [3]) but a different 
definition of nuclear shells Sj, where the first element is 
Si = J2j-rj<f ^j^ middle elements are Sj — for 
1 < rj < r, and the outer elements, rj > f, are Sj — Sj. 

In 2D and 3D the parameter -Bmax is proportional 
to the measurable maximum Overhauser field A = 
J2jLi^]^j^ ^//i is of the order of a few Tesla [l3 |. 



with the numerical factor tt ^ and (2Tr^ 



In ID, 



Bmax — NA/ (tt/i) is much larger than A, here N = 

J^'jLi'^SjAj/Ao. The parameter Bfluc = A^Js/N/n 
scales with the number of nuclei under the electron en- 
velop function in all dimensions. 

In terms of density of states the bare electron level 
acquires a finite smearing due to coupling to many de- 
grees of freedom of unpolarized nuclear spins. When the 
quantum dot is empty the nuclei at different lattice sites 



3 



are uncorrelated. After an electron, say with spin "up", 
populates the quantum dot, the state of the combined 
system I^^/q) — JJ^- --^ Ijl \]}-) is not an eigenstate of 
the Hamiltonian Eq. (|14p. where labels a subset of 
nuclear lattice sites and is the all spins down (includ- 
ing the central spin) state. We analyze the distribution 
of the eigenenergies Eq. ^ using a projected density of 
states = E{i^,k,}Pi{lj^k,})6{E-Ei{l,,k,})), 

where P{{lj,kj}) = 1 when (vJ/qI {Z^, fcj}) ^ and 
P{{lj,kj}) = when (*o|{^i,%}) = 0. Here the 
fc } runs over all subspaces and all eigenstates within 
each subspace. Note that for any shell with Sj > 1 the 
complete set of the eigenstates includes Ij with multiplic- 
ities greater than one Only one of each Ij is kept 
since these multiplicities do not change P {{lj,kj}). We 
calculate the overlaps matrix elements only in the leading 
l//i-B-order as the probability of measuring other eigen- 
states coming from subleading orders is at least as small 
as Aj/fiB. 

By representing the delta function as S (x) — 
J dAe""^/ (27r), the Fourier transform of ly (E) can be 
written as a product of sums over each nuclear spin shell 
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fc=Pj(l+sgnpj) 



(4) 



where pj — ('I'o|5'||^'o), \pj\ < Ij, are polarizations of 
the shells given by the state of the system |5'o)- 

Assuming that each shell is unpolarized pj <C Sj 
and Sj ^ 1, the sum within a shell can be calcu- 
lated as an integral, Ij (A) = /^f^ dfce~"'=('=+i) 



erf 
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1 + 25-, 



v^/2) -erf(v^/2) /{2Vix), 

X = XAj / (AjiB) , which is an oscillating function of 
A. Then the product of the oscillating functions 
can be approximated in the large-A^ limit by turn- 
ing it into an exponential of a sum of logarithms, 

UU ^0 W - h (A) exp (Er.r, >f log I, (A)) , and by 
expanding the exponent in A, j^r >f -^i (A) ~ 
Er.r,>A^^&S, - I {S,/2 + 5|/3) AA2/ (l^B) - (5|/24 + 
53/12 + 2Sf/45)X^Ay {iQfi^B^)]. 

In ID Ij (A) can not be calculated as an integral since 
the degeneracy of the hyperfine couplings is two but 
the explicit evaluation of the sum of only two terms 
within each shell and the small-A expansion yields a 
similar expression, J2f,r^yf log I j (X) « E^r,>f[log2~ 

iXA^^/ {AfiB) - }?Ay (a/24^B)^]. Strictly speaking, the 
small-A expansion is good when A ^ 16/ii?/>lp but the 
resulting Gaussian is also quite a good approximation for 
a large A since the original product of many oscillating 
functions is zero due to random phases of Ij (A) when 



A > 4:fj,B/Af, provided that the couplings Aj have a non 
regular distribution. 

By evaluating the inverse Fourier transform v (E) = 
J dXv (A) exp {—iEX) in the limit B 3> B^^i^y^ we obtain 



u{E)^ 



j:rj>f '-'3 



exp 



{E-Eo 



(5) 



where Eo = E™=iPj^j72 - iJ-B/2 is a shift of 
the bare electron level that depends on the mo- 
mentary state of the nuclei and a finite linewidth 



<^ = VS".,>f (^1/96 + 5^/48 54/90) ^ 

f^^Ruc/ (^y^B^ that is common for all unpolarized nu- 
clear states. In the intermediate regime Bfluc ^ -B <C 
Bmax Eq. (m is valid when E > SlAj/ (A^B). The 
contribution of the inner shells can be approximated as 
Ii (A) = Si when, due to the fast oscillating exponential, 
the main contribution to the inverse Fourier transform 

comes from A < 4:fj.B / (^SiAi^ . 

In ID, the Gaussian result agrees precisely with the 
spectroscopically measurable lineshape when B 3> Smax- 
As the degeneracy of hyperfine couplings is 2 for all 
shells, all projections |[13|] are the overlap of the singlet 
(or triplet) and two nuclear spin states which give I/V2 
and the calculation of the lineshape gives Eq. ([5]) . When 
the degeneracy is larger than 2 the two calculations are 
different. It is also worth noting that the state (0)) is 
an eigenstate of the model Eq. (1) with Sj = 1/2 in the 
high B-field B > B„,ax. 

Rediscretization of Eq. ^ recovers the average level 
spacing of the nuclear spin levels. From the definition of 
the density of states, d = l/u (E) is an energy range that 
contains only one state. But, as the prefactor in v {E) 
increases to infinity when more and more outer shells 
are taken into account, the level spacing becomes zero. 
On the other hand the coupling strengths of the outer 
shells become super-exponentially small which make the 
splitting of the inner shells' levels into sublevels due to the 
outer shells very narrow. Thus, by selecting an effective 
number of the significantly coupled nuclear shells rj < 
4ro, we find 

d(^) = rf(^o)exp[(S-^o)V^'] , (6) 

where d{EQ) = V^o"/ (^i nj.f<r,<4ro "^i) exponen- 
tially small, d{Eo) ~ exp (^-iV/^^ Thereby, d{Eo) 
is a tiny level spacing in the middle of the upper elec- 
tron Zeeman line but d {E) increases super-exponentially 
at a finite detuning E =/= Eq on a characteristic en- 
ergy scale a when B ^ -Bmax and SfAi/{AiJ,B) when 
Sfluc ^B <^B 

max- 

There is also a finite temperature smearing. To 
average the hyperfine shift Eq over all possible nu- 
clear spin configurations at a high temperature. 
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'^q{E) = T,{p^}HE - Eq), we use 
proach as in the calculation of 
the G aussian distri bution of levels with 

(To = 



the same ap- 
v (E) and obtain 
a width 



and an average level spacing do 



exp [{E - fiB/2f /al] Uj-.f <r,<4ro 2^.-/ iV^'^o)- 
This implies that if the nuclear spin state is not prepared 
in a specific way but is a thermal state, there are two 
energy scales in a projective measurement to narrow the 
nuclear spin bath |8l4l0l| in order to suppress fluctuations 
of the Overhauser field 0-01 • ^ measurement in the 
coarse resolution of do will select a single specific nuclear 
spin configuration suppressing only thermal fluctuations 
and a measurement in the fine resolution of d will project 
the system on an eigenstate within a given nuclear bath 
state. 

Using the eigenstates and the spectrum in Eqs. ([21 
|3]) one can evaluate the time-dependent density ma- 
trix of the electron with an unpolarized state of the 
nuclei, |vI/(0)) = {l + S+)U{J,^}It^m/V2 such that 
(\l/(0) |Jz|\l/(0)) = 0, as an initial condition. As a re- 
sult the diagonal matrix elements do not decay in time 
in the leading l/^B order, Ti = oo. When the degen- 
eracy of the hyperfine couplings is only 2 (fD case and 
B ^ -Bmax) the off-diagonal matrix elements have a slow 
Gaussian envelop with decay time T2 = l/cr on top of 
the fast electron spin Rabi oscillations with frequency 
liB. Note that one obtains the Gaussian decay assuming 
a phenomenological model of a quasistatic ensemble of 
nuclear magnetic fields Q- At a high temperature, av- 
eraging over different \^ (0)), one also obtains the Gaus- 
sian decay due to thermal fiuctuations with T2 = I/ctq 
[1^ which is much faster than l/cr. 

When the degeneracy of the hyperfine couplings is 
larger than 2 (2D and 3D cases and iJfluc <C -B <C i?inax 
in ID) we establish a bound on the shortest decay time 
assuming that all Clebsch-Gordon coefficients in the over- 
laps between the initial state (0)) and the eigenstates 
Eq. ^ are equal and neglecting degeneracies of Ij [l^l- 
This simplification gives a Gaussian decay with decay 
time T2 = l/cr. A more accurate calculation would give 
a spectroscopic lineshape, see discussion after Eq. ([SJ, 
which is narrower than the distribution of the eigenener- 
gies thus giving a longer decay time. 

The eigenenergies Eq. ^ are a good benchmark for 
numerical studies of Richardson equations [IJ. The spec- 
trum of the model Eq. (1) can be found at arbitrary 
field and for any quantum number n by solving a set of 
coupled non-linear equations Q, 



N 



A, 12 



1 - 



a.sE = E^ + hAj/2 - nB/2. At an infinitely 

large magnetic field solutions of these equations are sets 
of numbers E^, which are close either to —Aj/2 or ^iB. 



fiB 



2Ek 



— ^ El, — El- 



= 0, (7) 



At a finite magnetic field a 1/B-expansion of the Eqs. 
([15]) at these values of Ey recovers the 1 /i3-expansion in 
Eq. 121) and a 1/B-expansion of the Gaudin states 0j 
recovers Eq. ([3]). 

In conclusion we have diagonalized the central spin 
Hamiltonian in the high B-field limit. Projecting the 
eigenstates on an unpolarized state of the nuclear bath 
we have shown that the level spacing of the nuclear sub- 
levels, which is exponentially small in the middle of the 
bare electron level, becomes super-exponentially large 
with detuning away from the middle. This suggests to 
select states from the wings of the distribution when one 
attempts to eliminate the decohering effect of the nu- 
clei by projective measurement techniques. This theory 
is valid when the external _B-field is larger than typical 
Overhauser fields. 

We acknowledge support from the Swiss NF, NCCR 
Nanoscience Basel, DARPA, lARPA, and EPSRC. 



[I] I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 
323 (2004); R. Hanson, L. P. Kouwenhoven, J. R. Petta, 
S. Tarucha, and L. M. K. Vandersypen, Rev. Mod. Phys. 
79, 1217 (2007). 

[2] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

[3] A. V. Khaetskii, D. Loss, and L. Glazman, Phys. Rev. 

Lett. 88, 186802 (2002); Phys. Rev. B 67, 195329 (2003). 
[4] P. Simon, B. Braunecker, and D. Loss, Phys. Rev. B 77, 

045108 (2008). 

[5] W. A. Coish and D. Loss, Phys. Rev. B 70, 195340 (2004). 
[6] D. Klauser, W. A. Coish, and D. Loss, Phys. Rev. B 73, 
205302 (2006). 

[7] D. Stepanenko, G. Burkard, G. Giedke, and A. Imamoglu, 

Phys. Rev. Lett. 96, 136401 (2006). 
[8] A. Greilich, A. Shabaev, D. R. Yakovlev, Al. L. Efros, L 

A. Yugova, D. Renter, A. D. Wieck, and M. Bayer, Science 

317, 1896 (2007). 
[9] I. T. Vink, K. C. Nowack, F. H. L. Koppens, J. Danon, Y. 

V. Nazarov, and L. M. K. Vandersypen, Nature Physics 

5, 764 (2009). 

[10] H. Bluhm, S. Foletti, D. Mahalu, V. Umansky, and A. 
Yacoby, arXiv: 1003.4031 

[II] M. Gaudin, J. Phys. (Paris) 37, 1087 (1976). 
]12] See supplementary material. 

[13[ Representation of 2S-j > 2 spins 1/2 using a large an- 
gular momentum undercounts the number of the original 
states. The undercounted states are degenerate in energy 
and have to be accounted in a calculation of overlaps. 

[14[ D. Paget, G. Lampel, B. Sapoval, and V. L Safarov, Phys. 
Rev. B 15, 5780 (1977); M. Dobers, K. v. Klitzing, J. 
Schneider, G. Weimann, and K. Ploog, Phys. Rev. Lett. 
61, 1650 (1988). 

[15[ R. W. Richardson, Phys. Lett. 3, 227 (1963). 

[16[ J. Schliemann, A. V. Khaetskii, and D. Loss, Phys. Rev. 
B 66, 245303 (2002). 



5 



Supplementary materials 

Here we provide more details on the calculations and 
approximations of the main text. Figs. 1 and 2 illustrate 
construction of the shells in the Hamiltonian Eq. (1) of 
the main text. Section A contians a description of the 
zeroth order approximation to the eigenfunctions in the 
intermediate field regime. In section B we solve explicitly 
the Richardson equations in the large magnetic -B-field 
limit. 



Intermediate field regime Bn^c <C B -C Bmax 

The large magnetic field expansion has different con- 
ditions of applicability for the eigenenergies, 



N 



and the eigenstates, 



N 



m— 1 ^ 



(8) 



(9) 



in the subspaces of unpolarized nuclear spins kj 
Comparison of the first and the third terms in Eq. ([8]) set 
the limitation on external magnetic field as B^^^/ <C 1 
in all subspaces where 



Bru 



(10) 



The next order correction in Eq. ([5]) is small as B^^^/B'^ 
compared to the leading term. 

The second term of Eq. ^ is small when B ^ -Bfluc 
but the next order correction is not. The second order 
term in the flip-flop part of V contains a denominator 
which is a difference between unperturbed eigenenergies 
that belong to the same electron spin level. This sets the 
limitation on the magnetic field as 



[2^iB {A, - A, 



< 1 



(11) 



for all pairs i ^ j. Assuming that tq S> 1, the most 
restrictive condition comes from a pair of sequential shells 
at ri = 1 in ID and 2D as iJ/Smax ^ li where 



B„ 



2fi 



(12) 



In 3D the most restrictive condition comes from a pair 
of sequential shell at an intermediate radius rj = 
with 



8efi 



(13) 



where e is the base of the natural logarithm. In all cases 
^max is a much larger field than -Bfluc- 

The divergence of perturbation series for the eigen- 
states Eq. © when Bfluc <C -B ^ Bmax does not mean 
that a perturbation theory in 1/_B is inapplicable at all 
but only signals that the choice of the eigenfunctions, 

l*> - {S^f^'^^^U^=i {3+)"' 1^), is a poor zeroth or- 
der approximation to split the nuclear spin sublevels 
within one of the electron Zeeman levels. The correct 
approximation can be identified by merging the inner 
nuclear shells (in ID and 2D) with different couplings 
up to the radius f = rg \/ln (Bmax/ B) into a single shell 
with the same coupling Aq as the most restrictive limita- 
tion originates from the shells at the middle of the elec- 
tron envelope function. In 3D the correct approximation 
can be identified by merging the shells with intermedi- 
ate radii between f = t-q ^1 ± ^In {Bmax/ B)j /V2. The 

inner radius vanishes fast, when B = B^i^^/^, thus we 
neglect it and merge all of the inner shells up to the ra- 
dius f ^ ro (^1 + \/ln [Bn^sx/ B)^ /\/2 in 3D as well as in 
ID and 2D. 

Then we diagonalize H ^ V\ where V' = 
X]j r <f (^1 ~ ^i) ^0 • Sj, instead of the original model 
H from Eq. (1) of the main text by repeating the same 
calculation as for H and obtain Eqs. ([51 IH]) with a dif- 
ferent definition of nuclear shells Sj. where the first ele- 
ment is 5*1 — "Ylij r <f middle elements are Sj = 
for 1 < Tj < f, and the outer elements, rj > f, are 
Sj = Sj. Corrections to this result due to V are small 
for l^j — r| ^ 1. The specific form of the zeroth order 
eigenstates for rj ~ f can only be found numerically and 
we neglect this crossover region assuming a sharp transi- 
tion between the two types of eigenstates. 



A solution for the Richardson equations 



The eigenenergies of the central spin Hamiltonian, 
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H = fiBS^+J2^j^o-Sj 



(14) 



at arbitrary external magnetic field in each subspace with 
a given set of the quantum numbers n, L can be found 
by solving a set of Richardson equations [l| , 



N 

E 



2ljAj/2 



'■^^E^ + Aj/2 



2Ek 
En — Ek 



0, (15) 
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a 



N 



(16) 



All sets of E'l/S, which are the solution of the above set of 
equations, also uniquely define the Gaudin eigenfunctions 



\{EA) 




(17) 



that correspond to these eigenenergies. The unexcited 
state |JJ.) is all spins down, including the central spin, 
state. The normaliz ation f actor is determinant of annxn 
matrix 0|, z = V det M, which diagonal and the off- 
diagonal matrix elements are 

2lj {Aj/2y 



Mkk = 1 



E 



-[{Ek + A,l2y 



5Z 2i?i 
{Ek — Ep) 



Ml 



kk' — -2- (18) 

\Ek — Ek' ) 

This diagonlization procedure can be constructed in an 
easy way by solving a complementary bosonic model in- 
stead of Eq. to obtain the ansatz for the eigenstate, 
Eq. P?)l . Then, the Richardson equations emerge as the 
requirement for the states in Eq. (|17p to be the eigen- 
states of the model Eq. (IT4l) . In this way the exact form 
of Eq. (|15p can be found by using the spin commutation 
relations only. This approach was developed in [3] in the 
context of the BCS model and in [4] in the context of the 
Dicke model. 

When magnetic field is very large the last term in the 
Richardson equations, Eq. ([TSl) . can be neglected in lead- 
ing 1 / order therefore all of the roots are close to 
either —Ai,/2 or ^B. There is up to one root close to /ii? 
and, we consider here ID case only, there are up to two 
roots, kj < 2, close to ~Ai,/2. Collection of all possible 
sets of kj, such that n — fcj + (1 ± 1) /2 is equal to 

the number of excitations, explores the complete set of 
the eigenstates in a given subspace of the Hamiltonian. 
We will use —Ai,/2 and /ii? to label an eigenstate instead 
ofE,. 

Corrections to these limiting values of the roots at a 
finite B, E^ — —Ajjj2 -I- in and E^, — [iB -|- can 
be found from the 1/^iB expansion of the Richardson 
equations, Eq. p^ . 



5,= 




if one En=-An/2, 
if two Ev=-Ai,/2, 



if En = fiB. (19) 



The sum third case is over the remaining 

n ~ 1 roots —An/2. Note that the corrections to the 
roots —An/2 were obtained by linearizing the system of 
equations, Eq. (fT5| . and the correction to the root ^B 
was obtained by expanding Eq. ([15]) up to the second 
order in l/^B with the first order corrections to the roots 
—An/2 obtained from the linearized equations. 

The normalization factor z can also be expanded in a 
1/fiB series and we find that the product of the diagonal 
matrix elements of M gives in leading order 



'An/2, 



(2a)" riLi ir. if one En = ^iB, 



(20) 



where the product in the second case is over the remain- 
ing n—1 roots — An/2. Upper bound of the first sublead- 



ing correction to z can be estimated as rQyn\, assuming 



that a few off-diagonal matrix elements p8)) are of the 
order of Tq and they all contribute with the same sign to 
the determinant, which is smaller than the leading terms 
if ^B Aq \J n/e when number of the excitations is large 
^ 1 {Aq is the maximal coupling strength). We used 



the Stirling's formula to approximate 



for a large 



n. At high external magnetic field, B ^ i?fluc, the first 
subleading correction is small. 

Performing summation over En with the accuracy of 
Eq. (fT9l) we obtain the eigenenergies of the Hamiltonian 
Eq. (HH) at high magnetic field. 



E 



N 



2 ^ 

i=i 



Aj ( Ij + kj ) 



2L 



2 



2 ; 



A^jlB 



(21) 



where ± refers to the "up" and "down" electron Zeeman 
levels and kj labels the number of roots of Eq. ([T5t that 
are close to —An/2. Substituting En with the corrections 
from Eq. (fll?)) into Eq. (fT7|) we obtain the eigenfunctions 
that correspond to the eigenenergies E, 



N 



\{En})^{S+)^X{{S+Y' 14), 



(22) 



in leading l/fJ-B order. Corrections to this wave function 
contains admixture of states from the opposite electron 
level with single nuclear spin flip but each of them is 
proportional to a small factor Aj/ fiB. 

The l/_B-expansion of the Richardson equations in 
Eqs. (HIl [221) coincides with the result in Eqs. ([8l9)) 
obtained using the l/S-expansion of the main text. 
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Figure 2: Shells of constant couplings for a 2D square latt: 
by dashed lines select a subgroup of shells between r = 5 . 



;ice. The majority of the shells has degeneracy 8. Two circles marked 
and r = 6. 
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Figure 3: Distribution of shell spacings for a 2D square lattice, 100 < r < 101. The dashed line maxks the uniform spacing we 
use in the main text. 



